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Abstract 

We study the resonant x-ray scattering (RXS) spectra of CeBg in an antiferroquadrupole (AFQ) 
ordering phase, near the Ce L3 edge under the apphed magnetic field H || (1,1,0). On the basis 
of a locahzed electron model equipped with a mechanism that the RXS signal is brought about by 
the intra-atomic Coulomb interaction in Ce, we calculate the RXS spectra. The obtained spectra 
exhibit two contributions around the electric dipole {El) and quadrupole {E2) positions, and differ 
drastically when the orientation of H is reversed. The difference is brought about by the cross terms 
between the even-rank AFQ and magnetic-induced odd-rank contributions. At the El region, the 
relevant cross term is the one within the El process, while at the E2 region, they are ones within 
the E2 process and between the El and E2 processes. These findings capture the characteristic 
features the recent experimental data show, and provide a strong support and information of the 
field-induced multipole orderings. We also evaluate the RXS spectra near the L2 edge. Though the 
results show no E2 contribution, we find that the intensity around the El transition, which is as 
large as that at the L3 edge, can be detected experimentally. 

PACS numbers: 75.10.-b. 75.25.Dk. 78.70.Ck,78.20.Bh 
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I. INTRODUCTION 



Numerous and intensive research activities have been concentrated in the field of the 
/-electron systems where mutual interplay of charge, spin, and orbital degrees of freedom 
produces a rich variety of physical properties characterizing the strongly correlated electron 
systems. Among a large number of interests, the materialization of the multipolar ordering 
phase in these systems is one of the most fascinating topics to be addressed.-i^ 

Resonant x-ray scattering (RXS) is one of the most promising probes to observe exper- 
imental evidences of the multipole order with rank higher than two. A tensorial nature of 
the scattering amplitude helps our understanding of the experimental result." However, 
an analysis based merely on the symmetrical consideration of the tensorial nature does not 
answer the origin of the observed RXS signals. For instance, when the RXS signals were 
detected at the Mn K-edge in the orbital ordering phase of manganites.-i^ There was a 
controversy on what brought about the observed signals. In the electric dipole (El) tran- 
sition, the i^-edge signals reflect the anisotropic charge distribution of the Ap states, which 
form bands. Such an anisotropy may be given rise to by the distortion of the lattice or 
by the underlying ordering pattern of the 3d electrons through the intra-atomic Coulomb 
interaction.^ We call the latter as the "Coulomb mechanism" . Extensive investigations have 
revealed that the observed intensities were originated from the lattice distortion.—"— Since 
the 4p states of transition metals are rather extended in space, this result is quite rea- 
sonable. It is now recognized that the same mechanism is working on the RXS in other 
transition-metal compounds such as YTiOa and YVO3.— "— 

The situation is different for RXS near the L-edge in rare-earth compounds, since the 4/ 
states are so localized in space that the lattice distortion associated with multipole ordering 
is expected to be much smaller than that in transition- metal compounds. In CeBg, the 
antiferroquadrupolar (AFQ) ordering phase is inferred from various indirect observations 
such as macroscopic measurements, resonance methods, neutron scattering, and so on.— 
The most direct evidence of the AFQ order is provided by the RXS experiment by Nakao 
et al.'^'^ and Yakhou et al.^ who have succeeded in detecting the RXS signal at the Ce 
L3 edge at an AFQ Bragg spot G = (|,^,^)- Later, another experimental support was 
given by Tanaka et al. from the non-resonant x-ray Thomson scatteing (NRXTS) study 
which detected directly an evidence of the aspherical charge density. Note that no lattice 
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distortionis are observed in this material.— 

In our previous papers, we have calculated microscopically the RXS spectra on the basis 
of a localized electron picture which relies on the Coulomb mechanism, and have obtained the 
spectra in agreement with the experiments.— >^ In addition, invoking the same theoretical 



framework elaborated in Refs. 



27 and 



28 



we have calculated the NRXTS intensities on 
the AFQ Bragg spots. Our result has reproduced well the relative intensities of RXS and 
NRXTS and a Fano-dip like structure at G = (|, |, §),— in accordance with the Yakhou et 
aVs data.— 

It is predicted that the applied magnetic field induces the antiferro-octupole ordering in 
CeBg,— which could be detected by RXS and neutron scattering^^ measurements. Mat- 
sumura et al. have recently succeeded in detecting the octupole ordering induced by the ap- 
plied magnetic field via RXS.'^'^ By using different experimental settings examined before,— 
they observed the RXS spectra near the Ce L3 edge atG=(|,|,|) under the applied field 
along H II (1, 1,0) in the AFQ phase. Their data illustrate three notable features. First, in 
addition to main peak of the El transition around 5724 eV, the spectra show a small peak 
of the electric quadrupole (-E'2) transition around 5718 eV when the field is in the plus direc- 
tion (+H). This finding makes a remarkable contrast with the previous RXS data where the 
E2 peak is practically absent both experimentally and theoretically.— i^^iSiiS^ Second, when 
the orientation of the applied field is reversed, the spectral shape drastically changes. That 
is, the small peak around the E2 position becomes obscure when the field is in the minus 
direction (— H). Third, the peak intensity at the El peak varies in certain amount when the 
direction of H is reversed. Note that although the final feature was not emphasized in Ref. 
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, they confirmed that the difference actually exists, in particular, for the field strength 
approximately larger than 2 T.— 

These features were analyzed to originate from the cross terms between the contributions 
of primary AFQ (even rank) order parameters and those of field-induced (odd rank) order 
parameters. The purpose of this paper is to elucidate those observations by analyzing quan- 
titatively the spectra from a microscopic standpoint beyond qualitative one based merely 
on the symmetry consideration.^^ Developing our previous treatment,—""— we find that the 
field dependence comes from the interference effect between the even rank and odd rank sig- 
nals. Some of them come from the terms within the El process and within the E2 process, 
while others come from the terms between the El and E2 processes. Since the field-induced 



3 



octupole ordering, for example, gives rise to the third order rank signal, the information of 
octupole ordering could be extracted from the field dependence of the spectra. 

For numerical calculation, we adopt the same model and parameter settings to the pre- 
vious works as possible as we can to keep continuity of the research. We obtain the spectra 
reproducing semiquantitatively the three features the experiment had revealed. Note that 
the '/ast collision approximation\ which is widely used to analyze the spectra, is insufficient 
to discuss the field dependence, since it predicts no field dependence of the main peak in- 
tensity. Finally, we find that the RXS spectrum at the Ce L2 edge has enough intensity to 
be detected experimentally. 

This paper is organized as follows. In Sec. [Tll we introduce a theoretical framework 
to investigate the RXS intensity in the localized electron systems. In Sec. IIIH we briefiy 
summarize a model Hamiltonian which describes the initial state and used in our previous 
works, and explain the mean field solution of the Hamiltonian. Also the intermediate states 
of the scattering processes are presented. Numerical results of the calculated RXS spectra 
are shown in Sec. IIVI with comparisons with the experimental results. The last section is 
devoted to concluding remarks. 

II. THEORETICAL FRAMEWORK 
A. Scattering amplitude 

RXS is described as a second order optical process: photon with frequency u, wave 
number k, and polarization /i (= o" or tt) is diffracted by the sample into the state with 
the same frequency u, wave number k', and polarization fi' (= a' or tt'). The amplitude is 
approximated by a sum of the contributions from each Ce ion, which can be written as 

F(k, e; k', e':uj)=J2 /"(k, e; k', e'-.u), (2.1) 

n=l,2 

where /„ stands for the scattering amplitude of the En-En transition and electric 2"-th 
pole process is abbreviated to En. We omit the contributions from the terms like En-Em 
with n ^ m since we restrict our attention to centrosymmetric system. In this context, 
the En-En transition is simply called as the En transition hereafter. The En amplitude is 



4 



written as 

/„(k, e; k', e' : a;) oc ^ J] e-^^-r,^^^^- . ^ ^. ^ ^ . (3.2) 
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where Mn(j : k, e; k', e' : iS) represents the RXS amphtude at site j with position vector Vj. 
The number of Ce site is denoted as A^. Scattering vector G is defined by k'— k. Note that the 
above expressions aimed for absolute zero temperature are easily extended for a treatment 
of finite temperature (T) case by multiplying probability Pj^^ to M„(j : k, e; k', e' : iS) and 
summing over m. Here "Pj^m is proportional to the Boltzmann factor with Ej^rn being 

the energy of the m-th eigenstate at site j. For simplicity, we proceed the formulation for 
T = 0, though our numerical calculations in Sec. |IV]will be those for finite temperatures. 
The E\ amplitude is 
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(0|a:^,,|A)(A|V,,|0) 



fi,fi'=l A 

where |0) denotes the ground state with eigenenergy Eq, while |A) denotes the intermediate 
state with eigenenergy E\. The lifetime broadening width of the core hole is represented by 
r and it is fixed to be 1.5 eV in this work. The dipole operators described 
and Zj for fi = 1,2, and 3, respectively, in the coordinate system fixed to the crystal axes 
with the origin located at the center of site j. The E2 amplitude is 

/g2 \ 

M2(j : k,e;k',e : = y 2-] -eOvC^-e) 

(0|z,,|A)(A|v,|0) 
^nw-{E^-E,) + iT' 

where factors ^^(k ■ e') and g^'(k ■ e) with k = k'/|k'| and k = k/|k| are defined as a 
second-rank tensor. 







for fi = 


1, 






i(3A,5,- A-B) 


for fi = 


2, 




g/.(A,B) = < 


^{AyB,+A,By) 


for fi = 


3, 


(2.5) 






for fi = 


4, 






^ ~^AxBy + AyBx) 


for fi = 


5. 





Note that the quadrupole operator 5^ ,, is expressed as z^j = g^(rj , r^), and the subscripts 
fi = 1,2,3,4, and 5 for rank two quantity specify the Cartesian components — y'^,3z — 
r'^,yz,zx, and xy, respectively. 



In general, the evaluation of the scattering amplitude M„ tends to be a formidable task, 
since the intermediate states of the scattering process are difficult to calculate. However, 
there exist several cases in which the evaluation of the scattering amplitude becomes easy. 
For instance, by replacing the energy denominators in M„ with a single oscillator, the am- 
plitudes are reduced into compact forms.-*^ This treatment is called as ^fast collision ap- 
proximation'. Here, we adopt an another treatment, in which the Hamiltonian describing 
the intermediate states is assumed to preserve a spherical symmetry. For only the / states 
concerned, this assumption is justified when the crystal electric field (CEF) and the inter- 
site interaction are negligible compared with the intra-atomic Coulomb and the spin-orbit 
interactions in the intermediate states. Analyses based on this framework gave good results 
in many localized /-electron systems.—"— In the present case, this assumption seems appli- 
cable to the intermediate state of the E2 transition, while the applicability is not clear for 
the intermediate states of the El transition, since the 5d bands are involved. In this paper, 
assuming the same 5(i-D0S for the Cg and t2g symmetries, we preserve the spherical sym- 
metry in the intermediate states to calculate the RXS spectra. This assumption is justified 
later in our semi-quantitative analysis, since the RXS spectra are not sensitive to the shape 
and the fillingness of the 5d DOS. 

Within the present scheme, the RXS scattering amplitudes are expressed in neat forms 
suitable to discuss field dependence. For the El process, the scattering amplitude at a single 
site is 

2 2u+l 

M,{e-e' : u) = Y,a[^\u;) p(;)(e, 6')(0|zM|0), (2.6) 

u=0 11=1 

where zji''' is operator equivalence of mutipole moment of the component fi with rank u. 
For rank zero, zf^'^ = 1, and for rank one, zj^'^'s are Jx,Jy, and Jz with /i = 1,2, and 3, 

(2) 

respectively. For rank two, quadrupole operator is represented hj zli = g^(J,J). The 
ener gy p rofile of rank u contribution is denoted as a^^\uj), whose explicit form is found in 
Ref. 37. The geometrical factors are given as follows: for rank zero, P^(o)(e,e')= e' ■ e, for 
rank one, P^^^{e,e')= —i{e' x e)^, and for rank two, P^'^^{e,e')= g^(e',e). Note that we 
have omitted the subscript specifying the site and do so hereafter. For the E2 process, the 
scattering amphtude at a single site is 

1,2 4 2i/+l 

M,{k,e;k',e':u) = - ^ a^'^^M p(;)(e, e', k, k ) 

u=0 fi=l 
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X (0|4^)|0), (2.7) 
where P2'^^(e, e', k, k ) is the geometrical factor of the component fi with rank u. The def- 
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40 



initions of zjl^"^ with rank z/ higher than three and P2'^He, e', k, k ) are given in Ref. 
These formulae look similar to those derived in literatures, mainly employing the fast colli- 
sion approximation.-i^i^i^i^ Our treatment, however, is convenient when spectral analysis 
is needed since the energy profiles in Eqs. (12. 6p and (12.71) are correctly included. In the next 
section, we shall see Eq. (12. 7p is strictly applicable to describing the E2 process at the Ce 
L2.3 edges, while Eq. (12. 6p is approximately valid to express the El process. 



B. Intensity of the difference spectrum 

We investigate how the RXS intensity changes when the orientation of the applied field 
is reversed. For antiferro-type scattering vector G, the factor e"'^ ''^ appeared in Eq. (12. 2p 
gives +1 or —1 depending on the kind of sublattice the site j belonging to. By substituting 
Eqs. (O, ([22]), and ([22]) into Eq. ([21]), we obtain the total amplitude of RXS with H 
along a certain direction as 

P(k, e; k', e' : oc ViV Ua^^\uj)z[^^ + af ^(w)zf ^ 

+ia^^\oo)zi^^ + af{uj)zf + %af {uj)zf \ , (2.8) 

where the term ol2^(l>j)Z^^ is omitted since it is absorbed into (Xif^ {i>j)Z'^^ in the present 
case of CeBg.— Here, use has been made of a new quantity 

^ I Er=V -0(4^^) for V = even 

" 1 -^E';=iPrike,e'){zl;^^) for z/ = odd ' 

where the staggered moment is referred to as {zji''^). It is related to the sublattice moments as 
{zji'^)= {Olzj!^^ \0) A= —{0\zli'^\0)B where the subscripts A and B distinguish the sublattices. 
We emphasize that {zit^y^s are real quantities. Note that only staggered components of the 
multipole operators contribute to the amplitude with the antiferro-type G. 

The RXS intensity I{uj,H) is given by |P(k, e;k',e' : a;)p. When the direction of the 
applied field is reversed, z['^^ and Zg*^^ having odd u reverse their signs, while those having 
even u remain unchanged. Therefore, the amplitude for — H is expressed by the quantities 
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for +H. Then, the total intensities are expressed as 

I{uj, ±H) oc ai\uj)zf'^ + 

±i{a'l\uj)Z^^^ + ai\uj)zl^'^ + ai\uj)zf^] ' . (2.10) 
The difference spectrum is defined as 

The spectrum AI{uj) produces non-zero contribution when cross terms between the amph- 
tude with odd rank and that with even rank remains finite. By substituting fl2.1Up into Eq. 
f l2.1ip . we can classify such cross terms into three categories: 

A/(u;) = /\Ieie2 + ^Ieiei + ^Ie2E2, (2.12) 
Meie2 ^ 2Z^^^z!^hm[{a^^\uj)Ya^\uj)] 
+ 2Zf)z«Im[{a«(a;)}^af)(^)] 

+ 2Zf^zfhm[{af{uj)Yaf\uj)l (2.13) 

Meiei = 2Z^^^zfhm[{a^^\uj)Yaf\uj)l (2.14) 
AW = 2Z^^^zfhu,[{a^^\uj)ya^^\uj)] 

+ 2Zf^zfhm[{af{uj)yaf\uj)l (2.15) 

where Im[X] stands for imaginary part of X. 

Before showing numerical results of the RXS spectra, we comment on an outcome ex- 
pected from the fast collision approximation, in which the energy denominators in Eqs. 02.31) 
and (12. 4p are factored out of the summation over the intermediate states and replaced by a 
single oscillator. ^'^^ As a consequence, the energy profile an\oj) loses its v dependence, say 
an{oj). Then, AIeiei{'^) and AIe2E2{(^) become zero, and AIeie2{^) alone remains as 

AW(c^) = 2[Z«zf -Zf)(Z« + zf )] 

xlm[{ai{uj)ya2iuj)]. (2.16) 
This expression is what Matsumura et al. have used in their analysis.— 

III. INITIAL AND INTERMEDIATE STATES 

Cerium hexaboride is believed to show an AFQ ordering phase in the temperature range 
% < T < Tq with % = 2.3 and Tq = 3.3 K under no applied magnetic ^leld.^^^ 



It shows a simple cubic structure (CsCl-type, PrnSm) with lattice constant a being 4.14A. 
In CeBg, Ce is trivalent in the configuration forming a sextet term ^Fs. Under the cubic 
symmetry CEF, the sextet splits into a Tj doublet and a Fg quartet. The latter is the lowest 
energy state. Since the energy splitting between them is on the order of 530 K,^^ the Fg 
quartet alone is sufficient for investigating low-temperature phenomena. Within the Fg basis 
under the Oh symmetry, multipolar operators with rank one, two, and three are active. In 
the following, we introduce a model Hamiltonian defined in a subspace spanned by Fg basis. 



A. Model Hamiltonian 



In order to prepare the initial (or ground) state of the AFQ state in CeBg, we adopt 
the model Hamiltonian utilized in our previous works.—"— It is originally introduced by 
Ohkawa''^ and extended by Shiina et al.^^ It is derived on the basis of the RKKY interaction 
and a possibility of the anisotropic RKKY interaction is discarded for simplicity.— The 
Hamiltonian is 

5 2 



H 



/i=3 



(3.1) 



where operator equivalence of quadrupole moment is defined as = g^(J, J). The second 
line of Eq. f l3.ip describes the interactions between dipole and octupole moments, and the 
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and 
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The 



definitions of the symbols appeared in this line are found, e.g., in Refs. 
sum on runs over nearest neighbor Ce pairs. The last line in Eq. (13. ip stands for the 
Zeeman term with (/-factor being 6/7. Note that, by choosing parameter 6 positive, this 
Hamiltonian favors the AFQ order of three components belonging to the F5 basis {Oyz, Ozx, 
and Oxy)- The parameter 6 is fixed as 5 = 0.2. The coupling constant is chosen so as to 
reproduce the value of Tq in the absence of magnetic field. 
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FIG. 1: H-dependence of the staggered moments at T = 1.65 K. Bold solid and broken lines 
are (0-i,i,o) and (Jz), respectively. Thin solid and broken lines denote (t^) = {T^)/^/50 and 

\i'xyz 



B. Mean field solutions and the initial states 

We apply the mean field approximation to the Hamiltonian. Under the influence of the 
external field H along {h, k,i) direction, the mean field solution gives a ground state of the 



primary order parameter {Ohke) = {{hOyz + kOzx + (^Oxy)) /\/h? + k"^ + P. The quadrupole 
ordering temperature Tq in the zero field limit is given by 3.3 K for = 0.458 K. The field 
induces another ranks and/or another components of multipole order parameters depending 
on the direction of the applied field.— In the case of H || (1, 1, 0) adopted by Matsumura 
et al.,^^ three antiferro-type of the secondary order parameters, a dipole component Jz and 
octupole components T^yz and T" are induced. Note that the definitions of Txyz{= z^i^) 



and T"(= zf) are given in Ref. |38l The field dependences of these order parameters are 
shown in Fig. [T] for T = Tq/2 = 1.65 K. Since |H|-dependence of the order parameters, 
both primary and induced ones, become mild beyond |H| ~ IT, we will fix the field strength 
|H| = 2 T throughout the present work. This choice of the magnitude of H enables us to 
avoid complication concerning the problem of domain population, which is known to exist 
when the applied field is much smaller as seen in Fig. 2 of Ref. |33l Finally, note that though 
many more order parameters are induced in ferro-type alignments, we do not mention them 
since they have no contribution on the RXS intensity at the antiferro-type scattering vector 
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E1 process E2 process 

FIG. 2: Schematic diagrams of two RXS processes at the Ce L2 and L3 edges: (a) E\ process; 
(b) E2 process The 5d DOS is schematically exhibited in the upper part of the figure and arrows 
indicate spin of the bd states. 



now addressed. 



C. Intermediate states 



The intermediate states of RXS near the Ce ^2,3 absorption edges include excitations of 
an electron from core 2p states at a Ce site into the Sc? and 4/ states in the El and E2 
processes, respectively. Since the 5(i states form not levels but conduction bands, we need a 
model of the density of states (DOS) of them. In order to keep continuity from our previous 
work, we employ the same DOS used before.—^— The 5(i DOS, p^'^(x), is assumed to be 



0.008a; + 0.04, -5 < x < 0, 

O.Olx + 0.04, 0<x<8 (3.2) 
-0.0277X + 0.342, 8 < a; < 12.33, 



where x is measured in units of eV with x = corresponding to the Fermi level. Total number 
of the occupied electron per Ce site is set to be unity. We disregard the dependence on 
the bd states — y^, ?)Z^ — r^, yz^ zy, and xy. These settings are justified later in our semi- 
quantitative analysis, since the RXS spectra are not sensitive to the shape and the fillingness 
of the 5d DOS. Figure [2] shows a schematic view of the RXS processes and the shape of the 
5d DOS. 
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1. El process 



As explained above, the El transition at the Z/2,3 edges consists of that between the 
2p and Sc? states. We consider the resolvent [hw — iJj„t]~^, where Hint is the Hamiltonian 
spanned in the configuration involving one 4/ electron, one 2p core hole, and one excited 
electron in the 5(i band. First, we solve an eigenvalue problem considering the Coulomb 
interaction between the 4/ and 2p core hole as well as the spin-orbit interaction of them at 
the central site. Let the eigenvalue and the eigenstate be Ex and |A), respectively. In this 
work, the Slater integrals and the spin-orbit interaction parameters needed to evaluate the 
Coulomb and spin-orbit interactions are calculated for Ce'^+ atom within the Hartree-Fock 
(HF) approximation.— The obtained off-diagonal and diagonal values of the Slater integrals 
are multiplied by factors 0.8 and 0.25, respectively, taking the screening effect into account. 
Then, the presence of the bd electron is treated as a scattering problem. The problem is 
described by an inverse matrix problem symbolically summarized below: 

■ ^ dX,d'X' 



hw - Hint + il , 

[C'^hw + iT- Ex)y' SdXA'x' - VdX4'x' ^ , (3.3) 



-1 



where d = [rrid, Sd) specifies a state of the 5d electron. Matrix V stands for the Coulomb 
interactions between the 5d and 4/ electrons and between the 5d electron and the 2p core 
hole. The local Green function of the 5d electron G 5d{uj) is defined by 



G''\huj) = r —^^^ — de,7^0. (3.4) 



The RXS amplitude is calculated with rewriting Eq. (12. 3p as 



3 / I 

fl,^l'=l dx,d'X' ^ 



Hint / dXA'X' 



(d'A'|x^,,|0). (3.5) 



Detail of the resolvent treatment is found in Ref. 
Eq. ^ in Ref. 
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. Although the original derivation of 



371 does not expect inclusion of the hd band, the form of Eq. (12. 6p is 
still correct, since the bd DOS possesses a spherical nature. The dipole matrix element 
Adp={^d\r\2p) =J^ Rf^^{r)rR2p{r)r'^dr is included implicitly in Eq. (13. 5p where R2p{^) 
and i?gj(r) are the radial wave functions for the 2p and 5d states, respectively. Within the 
HF approximation, it is evaluated as Adp = 3.67 x lO^^^cm.— 
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scattered plioton 



FIG. 3: Schematic diagrams of the RXS geometry. The setting is adopted from Ref. [33 



2. E2 process 

The intermediate states in the E2 process can be constructed within the (2p)^(4/)^ con- 
figuration, disregarding the electrons preoccupied in the ground state. The Hamiltonian 
describing the intermediate states consists of the intra-atomic Coulomb and spin-orbit inter- 
actions in this configuration. The Slater integrals and the spin-orbit interaction parameters 
are calculated within the HF approximation.— The obtained off-diagonal and diagonal val- 
ues of the Slater integrals are multiplied by factors 0.8 and 0.25, respectively, taking the 
screening effect into account. The Hamiltonian matrix is numerically diagonalized by rep- 
resenting it in the (2p)^(4/)^ configuration. The RXS amplitude is calculated by inserting 
the eigenstates and eigenvalues into Eq. (12. 4p . Note that the scattering amplitude of this 
process is written by Eq. (12. 7p since the Hamiltonian preserves the spherical symmetry. 
In Eq. fl2.7p . the quadrupole matrix element Afp={Af\r'^\2p) =J^ R^^{r)r'^R2p{r)r'^dr is 
implicitly included where (r) denotes the radial wave function for the 4/ state. Within 
the HF approximation, it is evaluated as Afp = 5.69 x lO^^'^cm^.— 

IV. NUMERICAL RESULTS 



In the following, we shall report the numerical results of the RXS spectra. First, we 
clarify the setting of our RXS calculation. A schematic aspect is found in Fig. [3l in which 
the scattering vector G and the photon polarization are depicted. Contrary to our definition 
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FIG. 4: RXS spectra around the Ce L3 edges at G = (||^) with H || (TlO). Bold sohd and bold 
broken lines represent It^^u'{uj, +H) and I-^^ui{uj, — H), respectively. Inset: I{uj, ±H) around the E2 
peak. Thin solid and thin broken lines represent |a^^^(a;)p[zj^^]^ and \a^'^\uj) z[^^ -\- a'^\uj) Z^'^ , 
respectively. 



of G, some literatures, including the experimental works we analyze in the following,— >^ 
adopt the opposite sign, i.e., G^^P = k — k'. To avoid a confusion, when we mention G in 
this section, we mean k — k', while the actual calculations are carried out by using k' — k, 
because the formulae derived in our previous papers are the results of the latter definition. 



A. At G = (f |i) under H//(110) 
1. at the Ce L^-edge 

Matsumura et al. carried out the RXS experiment under the applied magnetic field 
H along (110) direction near Ce L3 absorption edge at G = (|||) in the vr-incident 
polarization.^^ They found the spectra showed strong enhancement of the intensity around 
the El and E2 regions. They also observed that the spectral shape underwent the significant 
change when the orientation of H was reversed. In particular, due to the cross term between 
the even rank and odd rank contributions, the spectra showed two-peak structure with peaks 
at the El and E2 positions when H was along (110), while they showed single-peak struc- 
ture with the E2 peak merged into the tail part of the El signal when H was along (110). 
We calculated the RXS spectra with the settings adjusted to those of Matsumura et aVs 
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(Fig. [3]). The results are shown in Fig. |H The origins of the energy are set and fixed such 
that the El and E2 peaks become around 5724 eV and 5718 eV, respectively. 

First, we concern the whole aspect of the spectral shapes. The calculated curves capture 
the above explained experimental feature well. We also confirm that the similar tendency 
is expected from the spectra in the a — n' channel (not shown) with the intensity of the 
El peak being roughly a fourth of that in the n — a' channel. Note that, the ratio of the 
intensity at the El peak to that at the E2 peak seems apparently quite different between 
the experiment and our result. For instance, for +H (|H| = 2 T), the ratio is about 2.45 
in the experiment^ and about 9.66 in our calculation. A remedy for this discrepancy is 
absorption correction. The experimental ratio is enhanced to about 10.5 after the correction 
is properly carried out.— So, our theoretical ratio gives fairly good values. 

Next, we analyze the ingredients of the intensity around the El and E2 positions. To 
this aim, we define average intensity as 

. /(^■+H)^/(..-H) ^ ^ 

At the El peak, by using Eq. (I2.10p . the average intensity is approximated as 

IM^) ^ \af\uj)nzi'^]\ (4.2) 

since another El contribution of |aj^^''((X')p[zj^'*]^ is two orders of magnitude smaller than 
that of |a;[^^(ci;)p[Z|;^-']^. At the E2 peak, the situation is not so simple. As shown in the 
inset of Fig. |U about half of the intensity around the E2 position is supplied by the tail 
part of the El contribution, |af''((X')| Other than that, an even-rank profile of the 

E2 transition (w) also has contribution by interfering with a\ (w) as illustrated by thin 
broken line. That is, 

/av(c^) ^ \c^?{oo)Z?^ + (a;)^f p. (4.3) 

Finally, we turn our attention to the difference spectrum AI{uj) displayed in Fig. O The 
sign around El peak and that around E2 peak are opposite to each other, which is also in 
accordance with the experiment.— Note that, as mentioned in Sec. [H although the AI{uj) 
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its existence is assured by more careful 



around the El region is not so prominent in Ref. 
measurement.— The presence of AI{u)) intensity at the El position, in itself, is ascribed to 
Eq. fl2.14p of AIeiei {Jz){Oiiq). It consists of the cross term between rank one and 
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FIG. 5: Difference spectra around tlie Ce L3 edges at G = (ff^) witli H || (110) in the vr — a' 
channel. Bold solid line represents the total A/(a;). Bold broken, bold dot-dashed, and thin solid 
lines are AIeiei{^^), ^lE2E2{i^), and AlEiE2ioj), respectively. 



rank two contributions arose both from the El transition, which is missing within the fast 
collision approximation as stated before. 

On the other hand, the origin of the intensity around the E2 peak is rather com- 
plicated. The whole shape is determined by A.Ieie2, while both AIeiei and A.Ie2E2 
have several quantitative contributions too. In the present setting, ^2^^ is two orders of 
magnitude larger than Zg^-* and the former is predominated by the contribution of T^yz- 
Then, AIeie2{^) and AIe2E2{<^) are approximated as 2Z[^^ Z^Hia[{a'^\uj)}*a^i\uj)] and 
2Z^^^ Z^^hui[{a'^\uj)}*a'^\uj)], respectively. Therefore, the entire spectral shape is well con- 
trolled by three energy profiles a^\uj), a^\uj), and a'^\oj) carrying by the AFQ component 
of Oj i Q and AFO component of T^yz- 

Here, we end this subsection with the explanations concerning the robustness of the 
present results. Among several assumptions we have employed in the present work, the 
choice of the So? DOS seems the most crucial one. The reason we have adopted the same 
DOS as used in our previous works is partly because to preserve the spherical symmetry 
in the intermediate states and partly to keep continuity of our research. Furthermore, we 
have confirmed that the characteristic features our data have shown are insensitive to the 
modification of the shape of the 5d DOS, fillingness of the 5d electron, and presence or 
absence of the DOS splitting between the Cg and states. For example, we have examined 
the semi-elliptic shape DOS and the uniform DOS. We also have changed the fillingness from 
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0.5 to 1.8 per Ce site and have introduced the DOS sphtting between the Cg and t2g states 
from —3 to +3 eV keeping the shape of the 5d DOS the same as that in Fig. [2J Although 
these modifications cause minor quantitative differences, the main features we stated in this 
section remain unchanged. 

2. at the Ce L2-edge 

Next, we also calculate the RXS spectra and the difference spectra expected from the 
AFQ phase in the vicinity of the Ce L2 edge where the El and E2 transitions are observed 
around 6167 eV and 6160 eV, respectively. The results are shown in Fig. O The intensities 
are nearly the same as or slightly stronger than those at the L3 edge (Fig. H]). We conclude 
that the RXS signal is experimentally detectable near the Ce L2 absorption edge. Then, 
we concern the spectral shapes. One notable feature is the E2 process has no practical 
contribution. This shows a striking difference compared to the RXS spectra detected near 
the Ce L2 edge in the antiferro-octupolar (AFO) phase of Cei-^^LaajBg where the signals 
from the E2 transition were distinctly observed.— i^"— The difference is attributed to that 
of the primary order parameter in these systems. 

B. At G = (lli) under H//(112) 

Nakao et al. detected the RXS signal of Ce L3 edge from the AFQ phase.— However, 
their result practically showed no E2 contribution contrary to the present case reported 
by Matsumura et al, in which the peak intensity at the E2 transition has distinguishable 
contribution from that at the El transition.— It is reasonable to interpret the difference is 
due to that of experimental conditions since the experiment by Nakao et al. was carried out 
at G = (^1^) under H/ /(112) and the polarization of the incident photon was a channel, 
all of which are different from the conditions chosen by Matsumura et al. 

In our previous works, we analyzed Nakao et aVs data and confirmed that the peak 
intensity at the E2 transition was negligible compared with the one at the El transition.— >2£ 
In these treatments, we set the peak position of E2 transition about ten eV lower than that 
of El transition. It turns out that the interval is too wide, since it is about six eV as 
seen from Matsumura et a/.'s data. If we set the E2 peak position at six eV lower than 
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FIG. 6: RXS spectra around the Ce L2 edges at G = (||^) with H || (110), respectively, (a) Sohd 
and broken Unes represent and I„fji{oj, — H) in the tt — a' channels. Inset: I(w,±H) 

around the E2 peak, (b) A/tto-' {(^) in the n — a' channel. 



that of El transition, both intensities may experience interference even under the Nakao 
et a/.'s experimental conditions. However, the calculated results (not shown), practically, 
have neither notable peak nor A/(ci;) intensities around the E2 position, which confirm our 
previous results survive after the shift of the E2 position. The absence of the distinct peak 
at the E2 position is merely the numerical reason. The tail part of the El contribution 
around the E2 position in Nakao et aVs case is nearly twice larger than that in Matsumura 
et aVs case. The former buries the E2 contribution, which is nearly the same in the latter 
case. 



V. CONCLUDING REMARKS 



We have theoretically investigated the RXS spectra observed in the AFQ phase of CeBg 
in the vicinity of the Ce L3 edge at G = under the applied field H along || (1, 1, 0) 

direction. The experimental data show small but clear contribution from the E2 process as 
well as that from the main El process. Also the interference between rank even and rank 
odd contributions from the El and/ or E2 signals are observed, which provides a great op- 
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portunity to obtain the information of the field-induced multipole orderings . To analyze the 
RXS spectra, we have employed the model on the basis of a localized electron picture, which 
is used in explaining several aspects of the RXS phenomena in CeBg and Cei_a;La2;B6.— ~— 
The model is combined with the intermediate states including the intra-atomic Coulomb 
and the spin-orbit interactions within the appropriate electron configurations, so that the 
RXS signal is brought about by the Coulomb interaction. 

The calculated spectra successfully capture the characteristic features the experimental 
data show: the ratio between the El and E2 intensities, the interference effect between the 
terms of rank even and odd when the direction of the magnetic field is reversed, and the 
signs of the difference spectrum AI{uj) around the El and E2 regions. 

When we focus on the detail of the spectra, the whole shape of the spectrum is roughly 
approximated by Jav = + a^\oo)Z^^\'^, i.e., contributions from the AFQ order 

parameter are dominant. Even around the E2 region, about half of the intensity is supplied 
by the tail part of the El signal. On the other hand, the difference spectrum AI{u)) is a 
direct consequence of the cross terms between the primary AFQ order and the magnetic- 
induced secondary order parameters with odd rank. A finite intensity of AI{u)) around the 
El peak, AIeiei{^), is observed by both the experiment and our calculation. This term is 
absent within the fast collision approximation. The AI{uj) around the E2 position consists 
of AIeie2{^) + ^Ie2E2{'^)- Comparing the A/(a;) with the fitting curve based on the fast 
collision approximation constructed by two Lorentzian curves in Ref. 31, we see that both 
give relatively similar outlooks. However, this is a coincidence. Our analysis have showed 
the main ingredients of the entire spectrum are four profiles a"2\'^)i 0^2^"* (i^)) and 
a"^\uj). We suggest that the difference of the spectral shape between our calculation and 
the outcome of the fast collision approximation may be observed when the orientation of 
the applied field is changed since the mixing ratios of four profiles are modified. A research 
toward this direction will be a future work. 

We also have found that the RXS signal is strong enough to be detected experimentally at 
the L2 edge, though the E2 peak is practically absent. We assert that the spectral analysis 
based on the microscopic calculation, which is beyond mere symmetrical consideration, is 
very useful to get deeper insights of the RXS phenomena. 

As mentioned in the preceding section, most of our results are robust semi-quantitatively 
when several modifications are introduced into the bd DOS. An exception is the fine structure 
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of AI{u). Utilizing the better (and/or realistic) 5d DOS may improve the detail of AI{u). 
There is a suggestion that the 5d DOS of the Cg and states are splitting and have the 
different shapes derived by a electronic structure calculation.^^ The splitting is expected 
at least in the high temperature region. We have checked, however, that such difference 
modifies little the overall shape of Ig,v{uj), while it affects in a subtle way that of AI{uj). If 
the latter spectrum would be measured with more precision, the spectral shape can be used 
to infer the form of the 5d DOS. 

A recent experimental report on magnetic spin resonance suggested that the AFQ-based 
model such as adopted in the present work met a difficulty in explaining its experimental 
data.— Even if a tiny amount of antiferromagnetic (AFM) moment is present, our results 
survive as long as the moment is small. A further test to our theory can be performed when 
the RXS signals are experimentally examined in the AFM phase. 
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